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Abstract 

The AdS/CFT correspondence defines a sector with universal strongly 
coupled dynamics in the field theory as the dual of pure gravity in AdS de- 
scribed by Einstein's equation with a negative cosmological constant. We 
explain here, from the field-theoretic viewpoint how the dynamics in this sec- 
tor gets determined by the expectation value of the energy-momentum tensor 
alone. We first show that the Boltzmann equation has very special solutions 
which could be functionally completely determined in terms of the energy- 
momentum tensor alone. We call these solutions conservative solutions. We 
indicate why conservative solutions should also exist when we refine this ki- 
netic description to go closer to the exact microscopic theory or even move 
away from the regime of weak coupling so that no kinetic description could 
be employed. We argue that these conservative solutions form the universal 
sector dual to pure gravity at strong coupling and large N. Based on this ob- 
servation, we propose a regularity condition on the energy-momentum tensor 
so that the dual solution in pure gravity has a smooth future horizon. We also 
study if irreversibility emerges only at long time scales of observation, unlike 
the case of the Boltzmann equation. 
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1 INTRODUCTION 



The AdS/CFT correspondence [ Q has opened up new vistas in understanding strongly 
coupled phenomena in four-dimensional conformal field theories. At strong cou- 
pling and large N, the dual classical theory of gravity admits a consistent truncation 
to pure gravity in asymptotically AdSs space and is described by Einstein's equa- 
tion with a negative cosmological constant. Any such solution of Einstein's equa- 
tion, which has a smooth future horizon, describes a dual state of the field theory 
at finite temperature and its dynamics. The final temperature of the horizon corre- 
sponds to the temperature of final equilibrium in the field theory. The dynamics of 
gravity is governed only by the five-dimensional Einstein's equation in this sector. 
The dynamics is therefore always universal, i.e. it is completely independent of 
the detailed matter content and couplings of the field theory. In fact any conformal 
gauge theory at strong 't Hooft coupling and large N, which has a gravity dual, 
must contain this universal sector [2 ]. This universal sector is practically important 
because it can describe, for instance, hydrodynamics of the dual conformal field 
theory (for a recent review, please see HI). 

However this universal sector, as described by gravity, covers a variety of phe- 
nomena well beyond the hydrodynamic regime. Such phenomena also include 
decoherence @. It would then be certainly useful to understand the universal sec- 
tor completely and also to decipher unambiguously the corresponding phenomena 
in the field theory side. This motivates our work. 

In the universal sector the dynamics of all gauge fields and higher form fields 
have been turned off on the gravity side. This means that in the dual field theory all 
chemical potentials vanish and therefore all conserved charge currents are absent. 
We strictly restrict ourselves to the case of zero global angular momentum in the 
field theory configurations, so that the final equilibrium configuration in gravity is 
always a static black brane. We will not explicitly mention this restriction in the 
rest of this paper; certainly most of the results here can be readily generalized to 
the case of nonzero global angular momentum. We also limit ourselves to the case 
of a flat boundary metric, so that the gauge theory is living in (the conformal class 
of) flat Minkowski space. The boundary topology then is R 4 and the static black 
brane is the unique equilibrium configuration in the absence of all other conserved 
charges. Metastable configurations like small black holes do not appear with this 
choice of boundary topology and boundary metric. 

The crucial aspect of the universal sector is that all the solutions in gravity are 
uniquely determined by the boundary metric and the boundary stress tensor. In 
particular, when the boundary metric is flat, we have good understanding [|6]-[H 
about the general features of the five-dimensional geometry [9 ]. It follows that all 
the states in the dual field theory constituting the universal sector are also uniquely 
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specified by the expectation value of the energy-momentum tensor alone. The 
dynamics of these states can then be determined by simply following the evolution 
of the energy-momentum tensor alone. This necessitates an understanding, from 
the field-theoretic point of view, as to how all observables and their time evolution 
could be functionally determined by the energy-momentum tensor. 

Here we will address this question; first in the regime of weak coupling, so 
that we can employ the quasiparticle description and also use kinetic theories, 
which are coarse-grained descriptions of microscopic laws. Specifically, we use 
the Boltzmann equation which has proven useful If 1011 1 1 H in determining the shear 
viscosity and higher order hydrodynamic transport coefficients and the relaxation 
time in weakly coupled gauge field theories. It has also been shown that an effec- 
tive Boltzmann equation can be used to study nonequilibrium phenomena in high 
temperature QCD and is equivalent to an exact perturbative treatment [10]. De- 
spite being a coarse-grained description, the Boltzmann equation retains the power 
to describe nonequilibrium phenomena far away from the hydrodynamic regime 
and at length scales and time scales shorter than the mean-free path and the relax- 
ation time respectively. However it is not applicable to phenomena at microscopic 
length and time scales. 

We prove that there exist very special solutions of the Boltzmann equation 
which are functionally determined by the energy-momentum tensor alone. We call 
such solutions "conservative solutions". These solutions, although very special, 
constitute phenomena far away from equilibrium and well beyond the hydrody- 
namic regime. The existence of conservative solutions can be conveniently proven 
for nonrelativistic monoatomic gases using some basic structural properties of the 
Boltzmann equation and can be easily extended to include relativistic and semi- 
classical corrections. We show that these solutions can be constructed even for 
multicomponent systems relevant for relativistic quantum gauge theories. 

It will thus be natural to make the assumption that the conservative solutions 
constitute the universal sector of strongly coupled gauge theories with gravity 
duals. This will explain why the states in the universal sector are determinable 
functionally by the energy-momentum tensor alone. This assumption, through the 
AdS/CFT correspondence, will have powerful consequences for gravity. The same 
condition on the energy-momentum tensor,required to make the state in the field 
theory a conservative solution, will now be required to make the dual solution 
in gravity have a smooth future horizon. In other words, the conservative con- 
dition on the energy-momentum tensor in field theory should now transform into 
the regularity condition in gravity. We will use this observation (assumption) to 
propose five equations which, in combination with the four conservation (or hy- 
drodynamic) equations, will provide a framework to determine the evolution of the 
energy-momentum tensor. Any energy-momentum tensor satisfying these equa- 
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tions will give us solutions in gravity with smooth future horizons. This framework 
has sufficient predictive power to determine nonequilibrium states in the universal 
sector beyond hydrodynamics in a systematic perturbative manner, given that the 
purely hydrodynamic energy-momentum tensor up to second order in the deriva- 
tive expansion (TJQJJEl is known. 

The plan of the paper is as follows. In Section 2 we outline the conservative 
solutions in the Boltzmann equation. We then state and investigate our proposal for 
the regularity condition on the energy-momentum tensor for pure gravity in AdSs 
in Section 3. Finally, in the Discussion we point out the various novel and open 
issues that have been raised in the course of development of our results and pro- 
posal here; particularly those which we may hope to understand in the near future. 
The proof of existence of the conservative solutions in the Boltzmann equation is 
slightly technical and elaborate; we present this proof in detail in the Appendix in 
a self-contained manner. 

2 THE CONSERVATIVE SOLUTIONS OF THE BOLTZ- 
MANN EQUATION 

The study of equilibrium and transport properties of dilute gases through the dy- 
namics of one -particle phase space distribution functions was pioneered by Maxwell 
lfl4l and further developed by Boltzmann lfl5l in the 19th century. The Boltz- 
mann equation provides a successful description of nonequilibrium phenomena in 
rarefied monoatomic gases. It is an equation for the evolution of the one-particle 
phase space distribution function. It can successfully describe nonequilibrium phe- 
nomena in rarefied gases, even at length scales between the microscopic molecular 
length scale and the mean-free path, and time scales between the time it takes 
to complete binary molecular collisions ifTTl and the average time between inter- 
molecular collisions. 

The Boltzmann equation is neither microscopic nor phenomenological, but a 
result of averaging the dynamics over microscopic length scales and time scales. 
Unlike phenomenological equations, it has no undetermined parameters and is 
completely fixed once the intermolecular force law is known. The structural de- 
tails of the molecules are however ignored and effectively they are taken to be 
pointlike particles. The hydrodynamic equations with all the transport coefficients 
can be determined from the Boltzmann equation. 

We start with a brief description of the conservative solutions of the Boltzmann 
equation for a system of pointlike classical nonrelativistic particles interacting via a 
central force. As mentioned in the Introduction, the proof of existence and unique- 
ness of such solutions is detailed in the Appendix in a self-contained manner. This 
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is followed by a discussion on how to generalize our construction of conservative 
solutions to the semiclassical and relativistic versions of the Boltzmann equation. 
Finally we show how our results apply to multicomponent systems relevant for 
relativistic gauge theories. These generalizations are straightforward and the dis- 
cussion on the nonrelativistic Boltzmann equation will be convenient for a first 
understanding of the conservative solutions. 

2.1 The conservative solutions in brief 

A generic solution of Boltzmann equation (1201 ) is characterized by infinite number 
of local variables. In general, these could be chosen to be the infinite local velocity 
moments (/ ( "^(x)'s) of the one-particle phase space distribution /(jc,£), given by 



where a - £ - w ; (x, t) with m,-(x, t) being the local average velocity. 

However the first ten velocity moments suffice to parametrize the energy-momentum 
tensor. The conservative solutions, which are determined by the energy-momentum 
tensor alone, are thus a very special class of solutions obtained when the initial 
value data satisfy certain constraints. 

Another special class of solutions to the Boltzmann equation is actually well 
known in the literature. These are the normal solutions, where the local hydrody- 
namic variables given by the first five velocity moments of / suffice to describe 
the solution even when it is far from equilibrium. Our conservative solutions are 
a generalization of these normal solutions. We review the normal solutions below 
before describing the conservative solutions. 

2.1.1 The hydrodynamic equations and normal solutions 

It is well-known that the first five velocity moments of the Boltzmann equation 
(1201 . obtained by multiplying with (1, £ 2 ) and integrating over give the hydro- 
dynamic equations as 




(1) 
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where the hydrodynamic variables (p, w,, p) are respectively the local density, com- 
ponents of local average molecular velocity and the local pressure of the gas de- 
fined in terms of the average root mean square kinetic energy. In terms of the 
velocity moments 



The local temperature is defined through the local equation of state, (RT = pip) 
ifloll . The shear-stress tensor ptj and the heat flow vector 5, (defined through 5,- = 
Sijkdjk) are related to the velocity moments by 



where c, - £ - It can be easily seen from the definition that PijSij - 0. 

The collision term /(/, f) (as defined in d2TI )) does not contribute when deriv- 
ing the hydrodynamic equations (0 from the Boltzmann equation. The first five 
velocity moments of J(f,f) are zero owing to particle number, momentum and 
energy conservation as proven in the Appendix. 

It must be emphasized that, in the hydrodynamic equations ©, the shear-stress 
tensor and the heat flow vector S are functionally independent of the hydro- 
dynamic variables. However there exist unique algebraic solutions to these and all 
the higher moments fr'. (x, t), which are functionals of the hydrodynamic vari- 
ables. These functional forms contain only spatial derivatives of the hydrodynamic 
variables and can be systematically expanded in the so-called derivative expansion 
discussed below. This leads to the construction of the normal or purely hydro- 
dynamic solutions of the Boltzmann equation, which we discuss below. For a 
generic solution of the Boltzmann equation, the higher moments of / will have 
explicit time-dependent parts which are functionally independent of the hydrody- 
namic variables. 

The normal solutions of the Boltzmann equation lfT8l - l20l have been exten- 
sively discussed in ETI . These solutions can be determined in terms of the five 
hydrodynamic variables (p,Ui,p) alone. They describe situations far away from 
equilibrium, such that observables which vanish at equilibrium do not vanish any- 
more but are functionally determined in terms of the hydrodynamic variables and 




(3) 




(4) 
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their spatial derivatives. The existence of such solutions follows from the exis- 
tence of unique algebraic solutions (as functionals of the hydrodynamic variables) 
to the equations of motion of the higher moments. The functional forms of the 
shear-stress tensor and the heat flow vector, for instance, are given by 

ij 2 r] 2 (D ( 1 \\ 

Pij - riaij +fii—(d- u)o-ij + p 2 — I — crij - 2 \o- ik cr k j - -6ijcr lm o- lm \\ 

+/?3^r UdjT - l$ipTj + UpdjT + djpdiT - pijd^Tj 

+p5^{d i Td j T-^6 u d l Td 1 T ] j + ... , (5) 

5; = xdiT + ... , 
with the convective derivative D/Dt = d/dt + u t di, and 

(Tij = djUj + djUi - ^djjd ■ u , 

p \5R 
rj - W) +... , X=—r} + ... , 

where tj and appearing as in the Navier-Stokes equation and the Fourier's law of 
heat conduction, are the shear viscosity and heat conductivity respectively. lf > is a 
specific function of the local thermodynamic variables determined by the collision 
kernel of the Boltzmann equation. The /3^s are pure numbers that can be determined 
from the Boltzmann equation. The time derivative in D/Dt can be converted to 
spatial derivatives using the hydrodynamic equations of motion; in fact, up to the 
orders shown above, we can assume that the Euler equation is valid and that the 
heat conduction is adiabatic. 

The functional forms can be expanded systematically in the derivative expan- 
sion, which counts the number of spatial derivatives present in the expansion. The 
expansion parameter is the ratio of the typical length scale of variation of the hydro- 
dynamic variables with the mean-free path. This is true for all the higher moments 
of /. The functional forms of pn and 5,- ©, when substituted into the hydrody- 
namic equations (O, give us systematic corrections to the Navier-Stokes equation 
and Fourier heat conduction respectively which can be expanded in the derivative 
expansion scheme. 

The hydrodynamic equations are now the only dynamical equations. The higher 
moments are given algebraically in terms of the hydrodynamic variables and their 
spatial derivatives. The phase space distribution function / is completely deter- 
mined by the hydrodynamic variables through its velocity moments. The hydro- 
dynamic equations thus form a closed system of equations and any solution of this 
system can be lifted to a unique solution of the full Boltzmann equation. 
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Stewart has shown [22] that such normal solutions exist even for the relativistic 
and semiclassical Boltzmann equations. 



2.1.2 The conservative solutions 

We are able to prove that a more general class of special solutions to the Boltzmann 
equation - which we call conservative solutions - exist. Here we outline these 
solutions, leaving the details of the proof to the Appendix. These solutions can 
be completely determined in terms of the energy-momentum tensor, analogous to 
the normal solutions being completely determined in terms of the hydrodynamic 
variables. The energy-momentum tensor (as shown later) can be parametrized by 
the first ten moments of / : 

• i) the five hydrodynamic variables (p, m,-, p), and 

• ii) the five components of the shear-stress tensor ptj in a comoving locally 
inertial frame. 

Importantly, for a generic conservative solution the shear-stress tensor is an inde- 
pendent variable unlike the case of normal solutions, where it is a functional of the 
hydrodynamic variables. 

These ten independent variables satisfy the following equations of motion 

+ —(pu r ) = , 
at ox r 

dui 1 d{p8 ir + p ir ) 
— — l- u r 1 = (J , 

dt dx r p dx r 

dp d 2 But \dS r 

^7 + -Z-(UrP) + ^(pS ir + Pir)^~ + ~ - , (6) 

dt dx r 3 dx r 3 dx r 

d Pij d dS ijr 1 dS r 
— h — — (U r pij) H Oij 

dt dx r dx r 3 dx r 

dui dui 1 du r 

+ T~Pir + T-Pjr - ^ijp rs — 

dx r dx r 3 dx s 



p,q=0;p>q;(p,q)±(2fl) 



where #|^' 9) are determined by the collision kernel in the Boltzmann equation, 
y and p indicate abstractly all the p and q indices of the moments f p) and p q \ 
respectively. 
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The above equations are now a closed system of equations, just like the hydro- 
dynamic equations were in case of the normal solutions. All the higher moments 
appearing in the above equations are given as functionals of the hydrodynamic 
variables and the stress tensor. These functional forms are unique and special alge- 
braic solutions of the higher moments of /. For instance, the heat flow vector can 
be determined from 

\5pRdT 3 Ldp ir „„ dT 2p ir dp \ 
£,• = —£- — + — —\2RT— + lRp ir — — ^- 1 + (7) 

25(2) dxi 25(2) \ dx r dx r p dx r ) w 

The functional forms of all the higher moments, as for the heat flow vector above, 
can be expanded systematically in two expansion parameters e and 6. The parame- 
ter e is the old derivative expansion parameter - the ratio of the typical length scale 
of spatial variation to the mean-free path. The new parameter 8 is an amplitude 
expansion parameter, defined as the ratio of the typical amplitude of the nonhydro- 
dynamic shear-stress tensor with the hydrostatic pressure in the final equilibrium. 

The closed system of ten equations (O are thus the only dynamical equations 
and any solution of this system can be lifted to a full solution of the Boltzmann 
equation through the unique functional forms of the higher moments. 

The normal solutions, being independent of the stress tensor, are clearly a spe- 
cial class of conservative solutions. There is another interesting class of conserva- 
tive solutions which are homogenous or invariant under spatial translations. The 
phase space distribution function / is a function of v only for these homogenous 
solutions and the hydrodynamic variables are constants both over space and time 
[this can be easily seen from ©]. The shear-stress tensor and consequently all 
the higher moments are functions of time alone. Such solutions have dynamics in 
velocity space only and describe relaxation processes. 

In a generic solution of the Boltzmann equation, the dynamics at short time 
scales is more like the homogenous class, where the initial one-particle distri- 
bution relaxes to a local equilibrium given by a local Maxwellian distribution 
parametrized by the local values of the hydrodynamic variables. At long time 
scales the dynamics is more like the normal solutions, where the system goes to 
global equilibrium hydrodynamically. Thus conservative solutions, despite being 
mathematically special, capture both relaxation and hydrodynamics which consti- 
tute generic nonequilibrium processes in a phenomenological manner. In other 
words, the dynamics of the energy-momentum tensor alone given by © captures 
both relaxation and hydrodynamics in a systematic fashion. 
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2.2 Relativistic and semiclassical corrections to conservative solutions 



The proof for existence of conservative solutions in the nonrelativistic classical 
Boltzmann equation can be readily generalized to its semiclassical and relativistic 
versions. This is because all the properties of the collision term / required for the 
proof of the existence of conservative solutions carry over to the semiclassical and 
relativistic versions as well. 

Let us consider the semiclassical version of the collision term which takes into 
account quantum statistics. This was first obtained by Uehling and Uhlenbeck ll23l 
to be 

J(f,g) = J J^,e)B(9,V)dfd€d6 , 

jtf.n = [f(x,og(x,f')ntmr)-f(x,og(x,aT(og(f')} , 

where the + sign applies for bosons, the - sign for fermions and s is the spin of 
the particles comprising the system. The final velocities £ and are determined 
by the velocties £ and before the binary molecular collision according to the 
intermolecular force law. Importantly, now J(f,f) vanishes if and only if / is 
the Bose-Einstein or the Fermi-Dirac distribution in velocity space for bosons and 
fermions respectively, instead of being the Maxwellian distribution ll24l . 

The proof for the existence of conservative solutions in the nonrelativistic clas- 
sical case does not require any explicit form of the collision kernel J. Only certain 
key properties suffice, as will be evident from the proof. We can pursue the same 
strategy with the semiclassical corrections as well 1125)1 . 

One has to employ the Sonine polynomials, which are generalizations of Her- 
mite polynomials, to find solutions of the required algebraic solutions of the higher 
moments as in li26l . The main objection could be that for the proof of existence of 
solutions, we use a theorem due to Hilbert which is explicitly stated for the non- 
relativistic classical /. However the details are exactly the same as that for con- 
structing the normal solutions. It has been seen that normal solutions can indeed 
be constructed in the semiclassical case [22J, so there ought to be no obstruction 
to the construction of conservative solutions also. Indeed, our proof shows that we 
can construct the conservative solutions given that the normal solutions exist. 

The generalization in the relativistic case again holds on similar grounds as 
above. It is more convenient to use a covariant description now. Normal solutions 
of the semiclassical relativistic Boltzmann equation have also been constructed 
ll22l . So there should be no obstruction in constructing conservative solutions as 
well. 
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In fact the same arguments could be used to state that any solution of the rela- 
tivistic semiclassical Boltzmann equation at sufficiently late times can be approxi- 
mated by an appropriate conservative solution, since the maximum speed of prop- 
agation of linearized modes increases monotonically ||27l as more and more higher 
moments are included. 

2.3 Multicomponent systems 

So far we have pretended as if our system is composed of only one component 
or particle. However gauge theories have many species of particles and internal 
degrees of freedom, hence we need to understand how to extend our results to 
multicomponent systems. 

Let us consider the example of N = 4 super Yang-Mills theory. In the weakly 
coupled description we need to deal with all the adjoint fermions and scalars along 
with the gauge bosons; all these particles form a SUSY multiplet. We note that in 
the universal sector all charge densities or currents corresponding to local (gauge) 
and global [SO(6)r] charges are absent. Similarly we should not have any multi- 
pole moments of local or global charge distributions, because in the gravity side we 
have pure gravity only. Therefore, most naturally we should have that all members 
of the N - 4 SUSY multiplet, distinguished by their spin, global charge and color, 
should be present in equal density at all points in phase space. So we are justified 
in our analysis in dealing with a single phase space distribution /. The Boltzmann 
equation we have considered above is obtained after summing over interactions in 
all possible spin, charge and color channels. 

The situation should be similar in any other conformal gauge theory. We can 
still treat the spin, color and charge as internal degrees of freedom owing to mass 
degeneracy even though the particles do not form a SUSY multiplet. In the absence 
of any chemical potential, there should be equipartition at all points in phase space 
over these internal degrees of freedom. This should be the most natural weak 
coupling extrapolation of the situation in the universal sector, dual to pure gravity, 
where gravity is blind to all the internal degrees of freedom of the particles. 

3 APPLICATIONS TO PURE GRAVITY IN ADS 

We will now argue that conservative solutions should exist even in the exact micro- 
scopic theory. In the exact microscopic theory, we do not make any approximation 
over the microscopic degrees of freedom and their dynamics unlike the Boltzmann 
equation, though an appropriate averaging over the environmental degrees of free- 
dom is required to get the final equilibrium configuration. 
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To begin with, consider the BBGKY heirarchy of equations [28 ] which de- 
scribes a hierarchy of coupled semiclassical nonrelativistic equations for the evo- 
lution of multiparticle phase space distributions. This hierarchy is useful for devel- 
oping kinetic theory of liquids. If the hierarchy is not truncated, then it is equivalent 
to the exact microscopic description. It has been shown that normal or purely hy- 
drodynamic solutions to the untruncated hierarchy exist. These solutions lead to 
the determination of viscosity of liquids which behave correctly as a function of 
density and temperature |29l . It is therefore likely that the conservative solutions 
also exist for this system which means they are likely to exist for the microscopic 
nonequilibrium theory of nonrelativistic classical systems constituted by pointlike 
particles. 

Experiments at the Relativistic Heavy Ion Collider (RHIC) suggest that the 
evolution of quark-gluon plasma (QGP) can be well approximated by hydrody- 
namic equations, very soon after its formation from the fireball QUI . Given that 
the perturbative nonequilibrium dynamics of hot QCD for temperatures greater 
than the microscopic scale A is equivalent to a relativistic semiclassical Boltzmann 
equation iflOl . we know perturbatively normal or purely hydrodynamic solutions 
exist for these microscopic theories. In fact any generic solution of the relativistic 
semiclassical Boltzmann equation can be approximated by an appropriate normal 
solution at a sufficiently late time. The quick approach to almost purely hydrody- 
namic behavior for the strongly coupled QGP in RHIC suggests that even nonper- 
turbatively a normal solution should exist which could approximate the late-time 
behavior for any generic nonequilibrium state. It is also true that not all transport 
coefficients of generic conformal higher derivative hydrodynamics can be defined 
through linear response theory. The plausible route of defining these higher or- 
der transport coefficients could be through the construction of normal solutions 
in nonequilibrium quantum field theories. Extremely fast relaxation dynamics in 
quark-gluon plasma similarly suggest that conservative solutions should capture 
generic nonequilibrium behavior. This is because in such systems the approach 
to the conservative regime, where the dynamics is given in terms of the energy- 
momentum tensor alone, should be faster than in weakly coupled systems, where 
even the corrections to Navier-Stokes hydrodynamics are hard to determine exper- 
imentally. 

If we accept that conservative solutions exist in the exact microscopic theory, 
it is only natural to identify the conservative solutions with the universal sector at 
large N and strong coupling in gauge theories with gravity duals. Such an iden- 
tification explains the dynamics in the universal sector being determined by the 
energy-momentum tensor alone. We emphasize, however, that the conservative 
solutions become universal only at strong coupling and large N. 

An appropriate AdS/CFT argument can also be provided for the existence of 
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conservative solutions for finite N and coupling. In such cases, we need to consider 
higher derivative corrections to Einstein's equation and the full ten - dimensional 
equations of motion. There is no guarantee of a consistent truncation to pure grav- 
ity anymore. However, we can use holographic renormalization with Kaluza- Klein 
reduction to five dimensions |f3T| to argue that we can readily extend the solutions 
in the universal sector, perturbatively in the string tension {~ 1 / V/l in appropriate 
units) and string coupling (whose N dependence is 1 IN). This can be done by turn- 
ing off the normalizable mode of the dilaton while keeping its non-normalizable 
mode constant, turning off the normalizable and non-normalizable modes of all 
other fields while keeping the boundary metric flat and perturbatively correcting 
the energy-momentum tensor to appropriate orders of the string tension and string 
coupling, so that the gravity solution still has a future horizon regular up to desired 
orders in the perturbation expansion. These solutions, again by construction, are 
determined by energy-momentum tensor alone. Our claim that the conservative 
solutions exist in the exact microscopic theory at any value of coupling and ,/V is 
therefore validated. 

The identification of conservative solutions with the universal sector at strong 
coupling and large N for conformal gauge theories with gravity duals allows us to 
create a framework for solutions of pure gravity in AdS with regular future hori- 
zons. We first study the parametrization of the boundary stress tensor which will 
allow us to make the connection with nonequilibrium physics. Then we will pro- 
ceed to give a framework for regular solutions, with the only assumption being the 
identification of conservative solutions with the universal sector at strong coupling 
and large ,/V. Finally we will make some connections with known results. 

3.1 The energy-momentum tensor and nonequilibrium physics 

A general parametrization of the energy-momentum tensor allows us to connect 
gravity to the nonequilibrium physics of conformal gauge theories. This parametriza- 
tion has been first applied in the AdS/CFT context in lfl2t . The energy-momentum 
tensor is first written as 

V = t(0)fiv + Kfiv , (9) 

where is the part of the energy-momentum tensor in local equilibrium. It can 
be parametrized in conformal theories by the hydrodynamic variables, the timelike 
velocity (w^) and the temperature (T), as 

t(0) MV = (nTfiAuftUy + rj MV ) , (10) 
and 7fy v is the nonequilibrium part of the energy-momentum tensor. 
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If we define the four velocity u M to be the local velocity of energy transport 
and the temperature T such that 3(7iT) 4 = w^w v f^ v is the local energy density, then 
in the local frame defined through u^, the energy-momentum tensor must receive 
nonequilibrium contributions in the purely spatial block orthogonal to the four ve- 
locity. This means 

M% y = . (11) 

The constraints in Einstein's equations impose the tracelessness and conserva- 
tion condition on the energy-momentum tensor so that 

<9% v = => & ({jtT)\4 Uii u v + V )) = -<9% v , 

Tr(t) = => Tr(n) = . (12) 

In the second equation above, the implication for the tracelessness for comes 
from the fact that the equilibrium energy-momentum tensor as given by (PTOl ) is by 
itself traceless. 

In the dual theory these conditions are satisfied automatically owing to the full 
50(4,2) conformal invariance. Note that the first of the equations above is just 
the forced Euler equation and can be thought of as the equation of motion for the 
hydrodynamic variables. 

We can reinterpret a class of known solutions of pure gravity in AdS as the du- 
als of the normal solutions in the exact microscopic theory at strong coupling and 
large N. These solutions are the "tubewise black-brane solutions" [13 ] which, in 
any radial tube ending in a patch at the boundary, are approximately boosted black 
brane solutions corresponding to local equilibrium and can be parametrized by the 
hydrodynamic variables corresponding to the patch at the boundary. These solu- 
tions can be constructed perturbatively in the derivative expansion. The expansion 
parameter, being the ratio of length and time scale of variation of the local hydro- 
dynamic parameters and the mean-free path in final equilibrium, simply counts the 
number of boundary derivatives. We can identify these solutions as duals of normal 
solutions because the nonequilibrium part of the energy-momentum tensor can 
be parametrized by the hydrodynamic variables and their derivatives alone. 

The complete parametrization of the purely hydrodynamic in any con- 
formal theory is known up to second order in the derivative expansion.In this 
parametrization, aside from the shear viscosity four higher order transport coef- 
ficients appear [12,13], which can be fixed by requiring the regularity of the future 
horizon giving us the tubewise black-brane solutions [13 ]. 

Let us define the projection tensor P^ v which projects on the spatial slice locally 
orthogonal to the velocity field, so that 
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The hydrodynamic shear strain rate cr^y is defined as 




(13) 



We also introduce the hydrodynamic vorticity tensor, 



(14) 



The purely hydrodynamic n^ v up to second order in the derivative expansion, 
for the tubewise black-brane solutions, with all nonvanishing transport coefficients 
fixed by regularity is 



Having identified the normal solutions in the universal sector with a class of 
solutions which could in principle be constructed up to any order in the derivative 
expansion, we will now naturally extend this observation to a framework which 
captures all regular solutions in certain expansion parameters. 

3.2 The complete framework 

In the hydrodynamic case we had four hydrodynamic variables, so the conservation 
of the energy-momentum tensor alone is sufficient to determine the evolution in the 
boundary. However in the generic case we need an independent equation of motion 
for 7fy v . 

The regularity condition must be an equation for the evolution of similar 
to the last equation of ©, This is because, as per our argument, the conservative 
solutions should be identified with the universal sector at large N and strong cou- 
pling. However Eq. © came from an underlying Boltzmann equation. At strong 
coupling, we have no kinetic equation to guide us because a valid quasiparticle 
description at strong coupling is not known even for N - 4 supersymmetric Yang- 
Mills theory. Moreover an entropy current cannot be probably constructed beyond 
the class of purely hydrodynamic solutions, hence we cannot use any formalism 
like the Israel-Stewart-Muller formalism [32] to guess an equation for n MV . This is 





+{\n2){nT) 2 {(T a ^ojJ + a^taj) + 0{d 3 u) 



(15) 
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because we should not expect a monotonic approach to equilibrium, as in the case 
of the Boltzmann equation, when we go to the exact microscopic description ||33ll . 

The safest strategy therefore, will be to use only the following basic inputs 
without resorting to guesswork. 

• The first input is that the equation for has to be conformally covariant 
because the dual gauge theory is conformal. 

• The second input is that the solutions in the purely hydrodynamic sector 
are known exactly up to second order in the derivative expansion and, be- 
ing identified with the normal solution, should be special cases of our com- 
plete framework. The equation for must therefore have ( fT5l) , the purely 
hydrodynamic energy-momentum tensor known up to second order in the 
derivative expansion, as a solution up to those orders. 

With only these inputs, we will be able to propose the equation for n MV only 
up to certain orders of expansion in both the hydrodynamic and nonhydrodynamic 
expansion parameters about the equilibrium state. However we should consider 
the most general equation for which satisfies the above criteria. The expansion 
parameters are again the derivative expansion parameter (as in the purely hydro- 
dynamic sector, but with the spatio-temporal variation of n MV taken into account 
additionally) and the amplitude expansion parameter, which is the ratio of a typical 
value of Tfyy divided by the pressure in final equilibrium. 

Our proposal then amounts to the following equation of motion for 7t^ v , whose 
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solutions should give all the regular solutions of pure gravity in AdS 5 : 



(1 - A 3 ) 



(u ■ 3)71^ + -7^ V (d ■ 



o - (^v 

2nT 



+ tt v V 



) (u ■ d)up 



(2 -In 2) 



7i^ v + 2(jtT) V 



-^ 3 (2-ln2)(ffrr 



(a • 3)^ + • u) 



+ ufo^j (u ■ d)up 



-A 4 (ln2)(nTy(cr atI 0J a v + (t^ojJ) 



1 



2X x (nTf \cT afl cT v a - -P^a^a; 
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2A 2 



(2 -In 2) 



+ - 



(2-ln2)(7rr) 3 

+0(^,ndn,d\,n 2 du,nd 2 u,d 2 n,d^uAdu){d 2 u)Aduf) , (16) 

where the 0(n 3 , ndn, ...) term indicates the corrections which lie beyond our inputs. 

The corrections can only include terms of the structures displayed or those with 
more derivatives or containing more powers of n^ v or both. We cannot say much 
about these corrections because for purely hydrodynamic solutions, they contribute 
to the energy-momentum tensor at the third derivative order only and the general 
structure of the hydrodynamic energy-momentum tensor at third order in deriva- 
tives is not known. The four /1,'s (i = 1,2,3,4) are pure numbers. Though we 
have not been able to specify their values, they are not free parameters. Once their 
values are fixed by regularity of the future horizon for certain configurations, they 
should give the complete framework for the whole class of regular solutions. 

As we have already mentioned, this equation of motion ([TBI for the shear- 
stress tensor has to be supplemented by the conservation of energy-momentum 
tensor in the form given in (fT2l so that we have nine equations for the nine vari- 
ables (including the hydrodynamic variables) parametrizing the general nonequi- 
librium energy-momentum tensor. The tracelessness of the energy-momentum ten- 
sor begets the tracelessness of n MV as in (PT2l) and this, as we have mentioned before, 
has led to the requirement that our equation of motion for n MV should be Weyl co- 
variant. 

This equation is thus a phenomenological framework for the universal sector 
as a whole up to certain orders in perturbation about the final equilibrium state. 
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This framework governs both hydrodynamic and nonhydrodynamic situations and 
goes much beyond linear perturbation theory. This is however, only valid within 
the universal sector. Beyond this sector we need many other inputs other than the 
boundary energy-momentum tensor to specify the boundary state or the solutions 
in gravity. 

3.3 Checks, comparisons and comments 

We will begin with a couple of comments. The first comment is that our Eq. (fT6l ) 
does not hold well at early times in the generic case. At early times the terms 
with time derivatives of various orders coming from the higher order corrections to 
our equation would become important. We will soon see the effect of such time- 
derivative terms in a simple example. We give an argument why such terms with 
time-derivatives must appear in the higher order corrections ll34ll . Any data at early 
times in the bulk, which will result in smooth behavior in the future, should get 
reflected in terms of an infinite set of variables in the boundary. The only way we 
can represent this in terms of the energy-momentum tensor alone is to include its 
higher order time derivatives in the initial data, so the equation for evolution of the 
energy-momentum tensor should contain higher order time derivatives. 

The second comment is that, in the particular case of boost-invariant flow, we 
have a better structural understanding of the hydrodynamic behavior at higher or- 
ders in the derivative expansion | 35| . We can, in principle, use our procedure to 
give a framework for general boost-invariant flows at late times. However we will 
leave this for future work. Moreover, the basic logic of our proposal is to use the 
purely hydrodynamic behavior as an input and then extend this to the complete 
framework. So our proposal and its extension at higher orders, by construction, 
reproduce the hydrodynamic sound and shear branches of the quasinormal modes. 

We now develop a straightforward strategy to check our proposal. We could 
look at simple nonhydrodynamic configurations first and construct the bulk so- 
lution perturbatively in the amplitude expansion parameter to determine some of 
the /i's. Once these have been determined, we can construct bulk solutions cor- 
responding to a combination of hydrodynamic and nonhydrodynamic behaviors 
perturbatively in both the amplitude and derivative expansion parameters and then 
check if the regularity fixes those /I's to the same values. 

The simplest nonhydrodynamic configurations are the analogs of homogenous 
conservative solutions of the Boltzmann equation we have mentioned before and 
which describe pure relaxation dynamics. Such configurations are homogenous in 
space, but time dependent and satisfy the conservation equation trivially. In such 
configurations the flow is at rest, so that if = (1,0,0,0) and the temperature T is 
also a spatiotemporal constant. The nonequilibrium part of the energy-momentum 
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tensor satisfies the following conditions 

(i) the time-time component ttqo and the time-space components 7To,- for i = 1, 2, 3 
vanish and 

(ii) the space-space components ttij for i, j - 1 , 2, 3 are dependent only on time. 
The above conditions on tt„ v result in the conservation equation being trivially 

satisfied. It follows from our proposal (fTBT ) that regularity in the bulk implies that 
itij satisfy the following equation of motion : 

dnij 2nT 1 - A { - A 2 ( 1 \ d 2 7r, 7 

(1 - h) ~dF + (2^2)^' - (2 _ ln2)(nT) 3 [™J ~ 3 Vwrfc] - 0(— ) . 

(17) 

If we look at the linearized solution, we have 



— ), T n = (1 -A 3 ) 
T n 2nl 



nij = Jlijexpi-—), t„-(1- A 3 ) , (18) 



where Mij is a spatiotemporally constant matrix such that tftifiij = 0. This implies 
that we have a nonhydrodynamic mode such that when the wave vector k vanishes, 
the frequency a> becomes purely imaginary and equals -it" 1 , i.e u = -h~ l as 
k — > 0. There is however, no such mode in the quasinormal spectrum of black 
branes |[36l . This makes us conclude that /I3 = 1, so that at the linearized level the 
only solution of (fT71) is = 0. 

However, at the nonlinear level we still have nonhydrodynamic solutions given 

by 

2nT I-A1-A2 / 1 c _ \ „A 2n ih 



(2 -In!)"- (2-ln2)(nTf\ nankj -f ijmmni i = 0( ^F ) ' (19) 

In fact, up to the orders explicitly given above, the equation is nondynamical and 
predicts that we should, at least perturbatively, have tensor hair on the black-brane 
solution in pure gravity in AdS. This gives us the simplest nontrivial test of our 
proposal and also a means of determining A\ + A 2 . 

In this connection, we also note that the possible second order in the time- 
derivative correction in ( fl9l implies that we need not have a monotonic approach 
to equilibrium as we should have in the presence of an entropy current. 

We end here with some comments on the issue of connecting our proposal 
with physics of quasinormal modes of the black brane. The linearized limit of 
the conservation equation, along with our proposed Eq. (fT6l ). supports at most 
three branches of linearized fluctuations. We have further argued that the third 
branch giving pure relaxation dynamics is not present. However, we know that 
the quasinormal modes have infinite branches of higher overtones other than the 
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hydrodynamic sound and shear branches. This naive comparison is somewhat mis- 
placed [ 37 ] because, as we know, nonlinearities do affect linearized propagation in 
quantum field theories. Since our equations are actually equivalent to the nonequi- 
librium field theory equation of motion of the state in the gauge theory, we must 
take into account nonlinearities of our equation in the propagation of the energy- 
momentum tensor before making any comparison. We leave this for future work. 

We would also like to mention here that it is only natural that the higher over- 
tones are more like resonances and are built out of the dynamics of the nine degrees 
of freedom of the conformal energy-momentum tensor, as it would have been sur- 
prising if infinite branches in the spectrum in the universal sector would have been 
blind to the microscopic details of the theory like the matter content and couplings. 
Our framework suggests that these infinite branches could be obtained from the 
nonlinear dynamics of the energy-momentum tensor. However we should exhibit 
caution here because although these nonhydrodynamic higher overtones of quasi- 
normal modes are indeed regular linear perturbations of the black brane, it is yet 
to be demonstrated that these can be developed into complete regular solutions of 
Einstein's equation nonlinearly. 

4 DISCUSSION 

We will mention some of the developments on which we would like to focus in the 
future. The first could be in the realm of early-time dynamics, especially in the un- 
derstanding of decoherence. This should be convenient because we can understand 
a lot by just considering the higher order corrections to the homogenous nonhy- 
drodynamic configurations which solve ( fl9l ). We have already observed the possi- 
bility of an oscillatory approach to equilibrium here. To uncover the physics, we 
need to compare with homogenous conservative solutions in quantum kinetic the- 
ories which can capture the physics of decoherence. We can test whether the same 
dynamics of the energy-momentum tensor in conservative solutions of quantum 
kinetic theories which captures decoherence, also gives rise to horizon formation 
in the bulk. 

A second important issue would be a better understanding of whether the hy- 
drodynamic limit of nonequilibrium dynamics always leads to generation of an 
entropy current. The tubewise black-brane solutions, which by our logic should 
constitute the normal or purely hydrodynamic solutions at large N and strong cou- 
pling, indeed demonstrate the existence of a family of entropy currents ll38l . How- 
ever, unlike for the normal solutions of the Boltzmann equation, or in the Israel- 
Stewart-Muller formalism, these entropy currents are not of the form sif, where s 
could be interpreted as the nonequilibrium entropy density. We hope to get a better 
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understanding of the physics of this entropy current by investigating the existence 
and form of the entropy currents in the normal solutions of untruncated BBGKY 
heirarchy which, as mentioned before, are solutions of exact microscopic dynam- 
ics. 

Finally, we have given a framework for general universal nonequilibrium be- 
havior in strongly coupled gauge theories with gravity duals. It would be inter- 
esting to see how much of this framework may apply to physics of quark-gluon 
plasma at the RHIC. 
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APPENDIX : PROOF OF EXISTENCE OF CONSERVA- 
TIVE SOLUTIONS 

We will now present the details of our proof for the existence and uniqueness of 
conservative solutions of the Boltzmann equation. In order to keep the proof rea- 
sonably self-contained, we give further details on the Boltzmann equation and how 
one can obtain the hydrodynamic equations seen earlier. We follow the notation 
of 112 111431 mostly for this part of the discussion. This will be followed by the proof 
in full detail. 



A.l A short description of the Boltzmann equation 

The Boltzmann equation for the one-particle phase space distribution f(x, £) for 
a gas of nonrelativistic monoatomic molecules of unit mass interacting through a 
central force is 

{ 1 r+Z- 1 r)K*,0 = J(fJ)M . (20) 
at ax 

where 

J(f, g) = f (/(x, f ')*(*, f ') - /(x, £>s(x, f)) m V)dfdedd , (21) 
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is the collision integral. (£, £*) are the velocities of the molecules before a binary 
collision and , ) are their corresponding velocities after the collision. The 
angular coordinates (9, e) are the coordinates related to the collision, and V = 
is the relative velocity with magnitude V. We assume that the collision takes place 
due to a central force acting between the molecules. 



Figure 1 illustrates the coordinates (6, e) used for describing the collision. The 
black dot in the center of the figure refers to the first molecule-the target molecule. 
The dotted line indicates the trajectory of the second molecule which we call the 
bullet molecule, with respect to the target molecule. The target molecule is placed 
at the center where its trajectory comes closest to that of the bullet molecule. We 
have drawn a sphere around the target molecule and n is the unit vector in the 
direction of the point of closest approach of the bullet molecule. The beginning of 
the trajectory asymptotes in the direction opposite V and the end of the trajectory 
asymptotes in the direction opposite V , which is the relative velocity t; - C after 
the collision. The co-ordinates (r, e) are polar co-ordinates in the plane orthogonal 
to the plane containing the trajectory of the bullet molecule and the target molecule 
as shown in the figure. The radial coordinate r is just the impact parameter as 




YJ 



v 



Figure 1 : The collision coordinates 
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shown in the figure. The angular coordinate 9 is the angle between n and the initial 
relative velocity V. Thus the unit vector n is determined by the angular coordinates 
9 and e. 

Solving Newton's second law for the given central force, we can determine r 
as a function of 9 and V, i.e. if the force is known we know r{9, V). The collision 
kernel B(9, V) is defined as 

dr(6, V) 

B(9,V) = Vr^-^ . (22) 

Finally the velocities of the target and bullet molecule are related to the initial 
velocities of the target and bullet molecule kinematically through 

i = 6-nKn-V) , 

£' = $+m(n-Y) , (23) 

so that V • n = V • n. 

This completes our description of the Boltzmann equation. When the molecules 
interact via an attractive or repulsive central force which is proportional to the fifth 
inverse power of the distance p between the molecules, we say the system is a gas 
of Maxwellian molecules. The simplification for Maxwellian molecules is that r is 
independent of 9 which can be seen from the fact that the trajectories of both the 
target and the bullet molecules lie on the circumference of a circle in the center of 
mass frame. As a consequence, B is also independent of 9. 

To proceed further we need to develop some notation. Let <p(g) be a function 
of We will call it a collision invariant if 

a>(£, e = m + 0(O - <t>(o - <p(?') = o . (24) 

Clearly there are five collision invariants - (l,£,-,£ 2 ) - which we will collectively 
denote as tj/ a . 

Let us also define, for convenience of notation, 

J(f.8) = J(f ' 8)+ 2 J(g ' f) ■ (25) 

Notation: We will use the following notation in the rest of this section. Let 
A (m) and be two tensors of rank m and n respectively, completely symmetric in 
all their indices. Then, 

• A (m) B ( ' l) will denote the symmetric product of the tensors so that it is com- 
pletely symmetric in all its m + n indices. 
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• A im) B (n) will denote a tensor of rank (m + n) where all indices except the i-th 
index in have been completely symmetrized. 

• We will use v as in A v to denote all the m indices in A 

• If Ajj and By are symmetric second rank tensors, then (AjjBy + + + ++) 
will denote the combination of all the six terms required to make the sum 
symmetric in its indices i, j, k and /. 

The above notations will hold even when A or B is a tensorial operator containing 
spatial derivatives. 

The hydrodynamic equations can be derived from the Boltzmann equation as 
follows. Using symmetry one can easily prove that 

fmj(f,g)df = \ f®({,f,f,{*')J(f,g)d{ . (26) 

Using (l24l) it is clear that if <p{%) is a collision invariant, that is <p(g) = tA a (£), 
then 



Mt)J(f,8)dt = . (27) 
A special case of the preceding result gives 

fffJS)Jif,f)d€ = ■ (28) 
The Boltzmann equation [on multiplying by if/ a (t;) and integrating] implies 

£ + 1; (/»•«) -° ' (29) 

where p a are the locally conserved quantities defined by 

Pa = \ tyafdi; . (30) 

These equations are equivalent to the hydrodynamic equations © once we 
make the identifications [po = P, Pi = pui (i = 1,2,3), p4 - (3p/2)]. 
The next few velocity moments, needed for later reference, are 

Pij = J(ciCj - RT8ij)fd^ , 

Sijk = J CiCjCkfdt; , (31) 

where c, = ^ - w,-. 
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A.2 The moment equations 

Multiplying both sides of the Boltzmann equation by higher polynomials of £ and 
integrating over we find the equations satisfied by the moments / (n) 's for n > 2 
to be 

(2) 

^ + ^( Ui fin) +f ^m + ^ f (n)_}_ f (n-V ) ^ = j(n ) ^ 

dt dxi V lJ J > i dxi' p J dxi 

where 

/"> = \ eB{ff[-ff x )deded^ , (33) 

is the n-th velocity moment of the collision kernel. 
It can be shown that 

DO 

/?>= J] tf$*<P>T)tf ) J? . 04) 

p,q=0;p>q 

with a particular simplification for B..' kl ' , which can be written as 

B i f.^ ) {p,T) = B^{p,T)5 ik 8 jl . (35) 

For Maxwellian molecules, there is yet another remarkable simplification that 
£(n,p,q)> s wt nonzero only if p + q - n. This happens essentially because the 
collision kernel B{6, V) in (|2TI) is independent of 6 in this case (for more details 
please see fl43l ). 

We will also denote fjj^ as Qiju and its explicit form will be useful. 

A. 3 Formal Proof of Existence of Conservative Solutions 

We now outline the proof that demonstrates existence of conservative solutions 
for the Boltzmann equation. The one-particle phase space distribution / will be 
functionally determined by the hydrodynamic variables and the shear-stress tensor 
(and their spatial derivatives). It must be emphasized that we proceed exactly along 
the same lines as used by Enskog in proving the existence of the normal (or purely 
hydrodynamic solutions) of the Boltzmann equation. 

The proof for the existence of normal solutions of the Boltzmann equation lfT8l - 
21] (first given in Enskog 's thesis) rests on the following theorem due to Hilbert 

EH. 

Theorem: Consider the following linear integral equation for g: 

J{f Q ,g) + J{g,f ) ) = f K , (36) 
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where J(fo,g) is defined through ([27]) and fo is a locally Maxwellian distribution. 
This equation has a solution if and only if the source term 7C is orthogonal to the 
collision invariants tjj a so that: 



I 



i/, a Kc% = , (37) 



provided the potential U(p) satisfies the condition that \ U{p)\ > 0(p~ n+l ) as p — > 
for n > 5. [That is, when the distance (p) between molecules vanishes, the abso- 
lute value of the potential should grow faster than (1/p) .] Further the solution is 
unique up to an additive linear combination of the tya 's- 

This theorem will be important in proving the existence of conservative solu- 
tions too, wherein we have to actually solve for the functional dependence on the 
hydrodynamic variables and the shear-stress tensor. For any conservative solution, 
we will just need to specify the initial data for the hydrodynamic variables and the 
shear-stress tensor. The only requirement will be that these initial data are analytic, 
because the functional dependence of / on the hydrodynamic variables and the 
shear-stress tensor will involve spatial derivatives of all orders. Clearly all normal 
solutions are conservative solutions, but not vice versa. 

The method of proof can be briefly outlined thus. We will extract a purely non- 
hydrodynamic part from the shear-stress tensor pn, and denote it as p ( " h \ This p^ 
will satisfy a simpler equation of motion which schematically reads (dp^ nh '/dt) = 
Z^Li c n (p( nh ') n , involving just a single time derivative [although the initial data for 
p^P can have any (analytic) spatial dependence]. The full shear-stress tensor py 
can be solved as a functional of the hydrodynamic variables and the p\" h \ One can 

functionally invert this to reinstate /?,-,• as the independent variable in place of p (n,i) 
and also determine the equation for p^. In the process we will see that there is 
an interesting class of nontrivial homogenous conservative solutions, where all the 
hydrodynamic variables are constants over space and time, while the shear-stress 
tensor is exactly Py , which is just a function of time. This class of solutions is 
thus purely nonhydrodynamic, representing equilibration in velocity space. 
The proof begins by writing the Boltzmann equation abstractly as 

£> = /(/,/) , (38) 

where 

df df 

£>=4+f-^~ , (39) 
ot ox 



and J(f, f) is as defined through (121b . 
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For a conservative solution, / is a functional of the nonhydrodynamic shear- 
stress tensor, p^ h \x, t) and the five hydrodynamic variables, namely w ; (x, t), p(x, t) 
and T(x, t). We expand / in two formal expansion parameters e and 8 such that 

oo oo 

/=2Z e " <5 '"^") • (40) 

n=0 m=0 

The physical meanings of the expansion parameters will soon be made precise. 
For the moment, if the reader so pleases, she can think of e as a hydrodynamic 
and 8 as a nonhydrodynamic expansion parameter. Following Enskog, we will also 
expand the time derivative in powers of e and 8 as : 

oo oo a(n,m) 

t^yy^iil . (4i) 
at ^-f^i at v 

n=\ m=0 

The above expansion of the time derivative might seem a little strange, but it will 
be necessary for us precisely for the same reason it was necessary for Enskog - the 
solutions of the equations of motion of hydrodynamic variables and pfj 1 ^ cannot be 
expanded analytically in e and 8, though their equations of motion could be through 
the subdivision of the partial time derivative. The proof will actually rely on the 
subdivision of the equations of motion just as in Enskog's purely hydrodynamic 
normal solutions and will not require the solutions to have analytic expansions [40]. 
This automatically results in a similar expansion for D, such that 

• For n > 1 and for all m 

D (n,m) s y y d (k ' l) An-k,m-l) + g/^m) (n > l;m = 0, 1,2, .,.) 

(42) 



k=l 1=0 



• For n = m = 0, 

D im - . (43) 

With the assumption that / is a functional of the hydrodynamic variables and 
Py h \ the time derivative acts on / schematically as 

dt ^5(v* P ) at i^a^ud dt 

y df W'T) y fl/ ^ 
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Above, schematically denotes k-th order spatial derivatives. Any time derivative 
acting on a hydrodynamic variable can be replaced by a functional of the hydrody- 
namic variables and the nonhydrodynamic shear-stress tensor by using the hydro- 
dynamic equations of motion. These functional forms have a systematic derivative 
expansion in terms of the number of spatial derivatives present and contain only 
spatial derivatives and no time derivatives. So the expansion of the time derivative 
in e is actually a derivative expansion, where the expansion parameter e is the ratio 
of the typical length scale of spatial variation of / and the mean-free path. This 
naturally "explains" (l42l . 

On the other hand, it will be seen that the time derivative of the nonhydro- 
dynamic shear-stress tensor can be replaced, using its equation of motion, by an 
infinite series of polynomials of the nonhydrodynamic shear-stress tensor. Thus 
the expansion of the time derivative in 8 as in (l43l : but we expand the solution of 
the equation of motion as an amplitude expansion with the expansion parameter 
8 identified as the ratio of the typical amplitude of the nonhydrodynamic shear- 
stress tensor with the pressure in final equilibrium. For the moment, these are just 
claims, to be borne out by an appropriate definition of the expansion of / and the 
time derivative. 

A.3.1 Subdivisions in terms of e and 8 

We outline here the expansion of the various quantities in the Boltzmann equation 
and the full Boltzmann equation itself in terms of the two expansion parameters 
£ and 8 and thereby arrive at various constraints that must be satisfied by these 
expansions. Our proof eventually will involve recursion while expanding the full 
Boltzmann equation in these expansion parameters. 

1. In close analogy with Enskog's original subdivision of /, we impose some 
further properties on the subdivision of /. 



• First we require, as in the case of normal solutions of Enskog and 
Chapman, that the hydrodynamic variables are unexpanded in e and 8 
and therefore are exactly the same as in the zeroth-order solution /(o,o> 
which will turn out to be locally Maxwellian. This is required because 
solutions of the hydrodynamic equations cannot be expanded analyt- 
ically in these expansion parameters though the hydrodynamic equa- 
tions themselves could be, as mentioned above. Therefore we should 
have 





(45) 
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where iff a are the collision invariants (l,£,f 2 ). It follows that 

Pa = J iff a f(p,0)dt; ■ (46) 

p a are the locally conserved quantities defined through d30l ). We may 
recall that these are just some combinations of the hydrodynamic vari- 
ables. 

We also require that the purely nonhydrodynamic part of the shear- 
stress tensor, pf^ has no expansion in e and 8, analogous to the hy- 
drodynamic variables. Being purely nonhydrodynamic, it determines 
f(0,m) for all m, i.e. the part of / which is zeroth order in e, but con- 
tains all orders of 5 in the conservative solutions. Since it vanishes at 
equilibrium, it is of first order in 8 and is given exactly by /(o.i)- More 
explicitly, for m > 2 and n = 0, we should have 



so that 



(c t Cj-RT6ij)f (f)/ „yde = (m>2) , (47) 



pff = f{ctCj-RT6ij)f m y^ . (48) 



2. The subdivision of the time derivative is defined next. Following Enskog, 
we impose on the time-derivative the condition that 



P" 







dt 

D {n ' m) ip a d% = 0; (n> l,m = 0,1,2,...) . (49) 



Using d45l ), the second equation above can be simplified to 



d (n,m) p a q 



dt 



+ iH {Mn-i,m)dt = (n > 1, m = 0, 1, 2, ..) . (50) 



Since the p a are a redefinition of the hydrodynamic variables, this above 
condition amounts to expanding the hydrodynamic equations in a particular 
way. From this expansion we know how each subdivision of the time deriva- 
tive acts on the (unexpanded) hydrodynamic variables. It is clear from (1441 
that if we now specify how the subdivisions of the time derivative act on 
p.. , we have defined the time derivative. Indeed, we have to solve for the 
action of the time-derivative because specifying this amounts to proving the 
existence of conservative solutions iflTTl . 
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3. The next thing is to note that the full shear-stress tensor p,j (just like any 
other higher moment) has an expansion in both e and 5. If we denote Spij = 
Pij - p{" h \ then for n > 1 



need not vanish. The expansion of 6p~ as a functional of the hydrody- 

namic variables and p~"!^ in e is the derivative expansion, with the power of 
e essentially counting the number of spatial derivatives (which act both on 
hydrodynamic variables and the nonhydrodynamic shear-stress tensor). The 
expansion in 6 is the "amplitude" expansion in terms of p^, which we may 
recall is first order in 8. 

4. On the basis of the above subdivisions one can now expand both sides of 
d38T ) and equate the terms of the same order on both sides. This enables 
us to write down the following set of equations that J(f, f) must satisfy for 
different values of (n, m). 

• For n - m - 0, substituting (@2]), (@3]> and (|40]> in ([38]> we get 




(51) 



^(/(o,o),/(o,o)) - 



(52) 



so that /(o,o) has to be a locally Maxwellian distribution. 
• Using the above fact, for n = and m > 1, we get 



J(J '(0,0), /(0,m)) + J(f(0,m),f(0,0)) 



,9(0,0) 



Bt 




(m > 1) 



(53) 



• Finally, for n > 1 and for all m 



J(f(0,0)>f(n,m)) + J(f(n,m)if (0,0)) ~ 




(54) 
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• The S („ >OT ) are given by, for (n + m > 2) 

n—\ m-1 

S (n,m) ~ 2_j ^(ffc 1 )' f(n-k,m-k)) 
k=l 1=1 

n-l m-1 

+ ^ J(f(k,Q),f(n-k,m)) + ^ J(f(0J), f(n,m-k)) (55) 
/t=l /=1 
n- 1 m- 1 

+ ^ J(f(k,m),f(n-k,0)) + ^ J(f(n,l),f(Q,m-k)) 
k=l 1=1 

+J(f(n,0),f(0,m)) + J(f(0,m),f(n,0))', (n + m) > 2 

and 

5(0,1) = 5(1,0) = . (56) 

A.3.2 A recursive proof 

With all of the above, we will now prove the existence and uniqueness of conser- 
vative solutions recursively. Recall that the key idea in this proof is to understand 
how the time derivative operator acts on the hydrodynamic variables and the 
nonhydrodynamic part of the shear-stress tensor, p., . We already know the ac- 

IJ 



tion of this operator on the hydrodynamic variables from Eqs.(|49l> and (I50I ). Now 
we will solve for the action of this operator on The action of the time deriva- 
tive, when expanded in e and 6, can be understood by analyzing the subdivisions 
of the Boltzmann equation given by Eqs. (l52l ). (l53l) and (l54l) . 

1. It is clear from d52l ) that at the zeroth-order in m and n, /( 0; o) is a locally 
Maxwellian distribution which is uniquely fixed by the choice of the five 
hydrodynamic variables (|46l ) and hence can uniquely be specified as 

/(0,0) = ^^ ex p(-^W) • ( 57 ) 

(2nRTy2 \ 2RTJ 



2. Next let us consider (1531 . The usual trick here is to rewrite /(o, m ) as /(o,o)fyo,m)- 
The advantage is that since /(o,o) contains hydrodynamic variables only, 

g(0,m) 

— / ( o,o) = . (58) 
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Therefore (l53l) can be rewritten as 



5(0,0) 

= /(0,0) 2_j ~7~"(0,m-0 - 2^ ^C/(0,0)A(0,/)./(0,0)A(0,m-o)5 ( m ^ 2 )- 
Z=l /=1 

Now we analyze d59l for m = 1 and m — 2. 



J(f(0,0), f(P,0)h(0, m )) + J(f(0,0)h(0,m), /(0,0)) ~ /(0,0) —^-h(0,m) (59) 

3(0,/) m - 1 
/ 1 ~ 



m=l: 



For m-\, (1591 reduces to 



J(f(0,0), 7(0,0)^(0,1)) + ^(/(0,0)^(0,l)»/(0,0)) 

3(0,0) 



while it follows from (f48T> that 



= /(0,0)-^%i) , (60) 



%d = 2! p/?r foCj " * rtfy) " (61) 



These two equations imply that 

3(0,0) 



(H p ™ = BF>(p,T) Pl j , (62) 



where B (1) has been defined in ([35]) [42J. 
• m=2: 

At the second order, (l59l implies 



^(/(0,0)>/(0,0)fyo,2)) + A/(0,0)fyo,2),/(0,0)) 

3(0,0) 3(0,1) 
= f(P,0)—Q^-h(p,2) + /(0,0)-^y- %,i) - /(/(ao)%4)»/(0,oAo,i)) ■ 

(63) 

We then need to solve for two things, ft(o,2) and the operator fS^'^/Sf). 
To do this we first write /i(o,2) as 
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1 s 



1 



*(0,2) 



Q 



(0,2) 

m 



3! p(RT) 2 



(cicpk - RT(aS jk + +)) 



(64) 



4\(RT) 2 



ijkl 



p(RT) 



-5 M + + + ++ 



(SijSti + +) 



x[acjC k ci - RT(c iCj dki + + + ++) + (RT) 2 (6ijS k i + +)] . 

The idea behind guessing this form is to expand /z(o,2) in two higher 
order Hermite polynomials of c's and reexpressing the Hermite coef- 
ficients through the ordinary moments. This method of expansion is 
due to Grad [43 ]. For the moment we can just take it as the most gen- 
eral possible form of ft(o,2)» since if higher Hermite polynomials are 
included here, the coefficients would have vanished. It also turns out 
that sfj'^ vanishes. Similarly all the other higher odd moments vanish, 
so far as their purely nonhydro 'dynamic parts ( or expansion in m for 
n = 0) is concerned. Obviously this does not mean that these higher 
odd moments have no dependence on p^K For n > there is indeed a 
nonvanishing expansion in m for these moments. We can now compare 
the coefficients of Hermite polynomials on both sides of our Eq. (l63l . 
For Maxwellian molecules (thus determining the form of J) we have 



d ' (nh) _ R (2,2,2). T , (nft) (nft) R (2,4,0). T , n (0,2) „~ 
~QfPij ~ B ijklmn { fi' T >Pkl Pmn + B ijklmn^^'Qklmn ' (65) 

d (Q ' Q) „(Q,2) _ R (4,4,0) , _ T ^(0,2) , R (4,2,2) , T ,(nh) (nh) 
Q t ^ijkl ~ D ijkImnpqW> 1 >^mnpq t ijkhnn pq ' 1 Wmn Ppq 

Since we know the action of (d (0 '°V<9f) on p^ and the hydrodynamic 

variables, we can solve for Qfj^ as a functional of p "® and the hydro- 
dynamic variables; the solution turns out to be 

Q-klmn ~ XklmnpqrsPpq Prs j (66) 



where Xkimnpqrs satisfies the equation 1441 

9R(2) y _ R (4,4,0) y , d(4,2,2) (frj , 

A klmnpqrs ~ n kImnijtu Ai j !u P1 rs + D klmnpqrs ' ^ 1 ' 

This in turn provides the solution for the operator (d^'^/dtp 
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d (0,1) (nh) _ R (2,2,2) (nh) (nh) R (2,4,0) y (nh) (nh) , fi8 \ 

dt P 'j ~ 'jklmnPkl Pmn + a ij k i mn A klmnpqrsP 'pq Prs ■ 

The equation above shows the action of the operator on p {nh) ; we al- 
ready know how it acts on the hydrodynamic variables. This implies 
that we have solved for this operator at this order. Note that the solution 
for the operator corroborates the intuitive understanding that this oper- 
ator is the next order in amplitude expansion. Another important point 
is that the solution of the operator is not independent of the solution 
for and is just given by the logic of our expansion once QjjjJ has 

been solved as a functional of pf^. This feature is the same for all the 
higher terms in the expansion of the time derivative operator as well. 
For non-Maxwellian molecules things are a bit complicated because 
the equation for Qiju in (|65T > also contains a term linear in p ( "j h \ so that 
now 



d^°' l) (nh) c D (2) (nh) , D (2,2,2) (nh) (nh) 

—p\j = SB ( >p\. + B)j£y u } p K J 

R (2,4,0) Y n (nh)(nh) (f.Q\ 

ijklmn \<dmnpqrsPpq Prs ■ v u "/ 

However this feature also appears in the usual derivative expansion (the 
expansion in e) of the time-derivative. Despite appearance, SB^p^ 
is a small quantity as (SB^ '/B™) is a pure number which is smaller 
than unity (for a proof of this and also for the statement of convergence 
of such corrections in the context of normal solutions, please see lfT9l 
20]). This result can be translated here, as the normal solutions are 
just special cases of our conservative solutions and at a sufficiently late 
time our solutions will be just appropriate normal solutions P31 . 
This is indeed remarkable considering that we have no parametric sup- 
pression here. Formally however, aside from the convergence problem, 
there is no obstruction because 6 is just a formal parameter and is only 
intuitively connected to the amplitude expansion. 

Higher m: 

We can proceed in the same way to the next order in m when n is 
zero. At every stage we have to deal with /(o, m ) which we may write 
as f(p t o)h(p t m) and further expand /j(o, m ) in a series containing up to m-th 
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order Hermite polynomial in c's. We have to solve for the coefficients 
of these Hermite polynomials, which depend on x only, and this leads 
to the definition of the ra-fh subdivision of the time derivative operator 
in the 6 expansion when the e expansion is at the zeroth order. The 
equation for evolution of p (nh \x,t) thus finally involves only a single 

IJ 

time derivative which we have expanded in 5. This is highly nonlinear, 

(nh) 

involving an infinite series of p%". The presence of just a single time 

derivative in the equation of motion for p^" h \x, t) makes it essentially 
an ordinary differential equation in one variable and so for any initial 
data existence and uniqueness of solution is guaranteed. 
We note that we can consistently truncate our solution at n = so that 
there is no expansion of / in e, provided all the hydrodynamic variables 
are constants over both space and time and is constant over space 
but a function of time. This gives us the simplest class of conserva- 
tive solutions which is homogenous in space; the Boltzmann equation 
becomes equivalent to an ordinary differential equation involving a sin- 
gle time derivative for pij. Physically this solution corresponds to the 
most general conservative solution which is homogenous in space, but 
generically far away from equilibrium in the velocity space (so that the 
velocity distribution is far from being Maxwellian). 

3. The next task is to prove the existence of solutions for the recursive series of 
equations in (l54l ). To see if solutions will exist we need to employ Hilbert's 
theorem. S(„, m ) contains either pairs of the form J(f(p tq ), f(r, s )) + J(f(p,q), f(r,s)) 
or just /(/(/,/),/(/,/))■ So when the collision invariants are integrated with 
S(n,m), as in J if/ a S ( n ,m)d^, the integrals vanish as a consequence of (l27l) . 
Therefore, the existence of the solution to f{n, m ) follows from Hilbert's theo- 
rem as a consequence of ( |49l ). The solution is unique because the condition 
(l45l fixes the arbitrariness of the dependence of f( n ,m) on the collision in- 
variants \f/ a . The details for n > 1, are thus, exactly the same as in the case 
of normal solutions. The action of (<3 (,V "V<9?) on p..^ is also determined as 

' J 

soon as the functional dependence of Spij and the relevant higher moments 
on and the hydrodynamic variables are determined. 

The explicit calculations become extremely complex even when, say n = 
2,m = or n = l,m = 1. We give some explicit results for the first few 
terms in the expansion for Sptj as 
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~ R (2,2,2) 

L P a ijUmn (nh) ( dUm dUn_ _ 1_ r 
(B (2) ) 2 % ^ + 5Xm 3^"^/ • 

It is clear that the terms in the expansion involve spatial derivatives of both 
the hydrodynamic variables and p~ nh) . From the expression for 6p^ 1 ' ^ one 
can determine the shear viscosity r/ which is of course the same as in the 
purely hydrodynamic normal solutions, so that 

n*J^(p,T) . (7i) 



p du m du n 2 du 



We also give some terms in the expansion for the heat flow vector 

\5pRdT_ 

25(2) dxi 



5 a,o) = I5g5r (?2) 



5 (U) = ^\ 2RT d Ii^ +1Rp inH)^_ 2 _Pu_dP_\ 

dx r lr dx r p 8x r 



It is clear that the heat conductivity x is a l so the same as in purely hydrody- 
namic normal solutions so that 

\5R p \5R 
XK ~B^ (P,T)K ~ 11 ■ (73) 

Corrections to the above relation appear in the higher order for non-Maxwellian 
molecules but again these are the same as in the case of normal solutions. 

This completes our proof for the existence of conservative solutions for the 
nonrelativistic Boltzmann equation. As mentioned before, we can now reinstate 
Pij as the independent variable. Our independent variables satisfy the following 
equations of motion 



36 








dui dui 1 d(p5 ir + pir) 












(74) 



W Ptr + d7 Pjr " 3 SijPT 'dT 




TfL,p,q), T , Ap) At 
B iivp (P» y )/v /p 



p,9=0,P>9;(P.9)*(2,0)) 



The first three equations are just the hydrodynamic equations, while the equa- 
tion for pij can be obtained from (l32l . 

The crucial point of this proof is that we have now solved for all higher mo- 
ments / v s for n > 3 (which includes, of course, S ijk and thus 5,) as functionals 
of our ten variables (p, Ui,p,pij) with T = p/(Rp). Any solution of these ten equa- 
tions of motion can be uniquely lifted to a full solution of the Boltzmann equation 
as all the higher moments are dependent on these ten variables through a unique 
functional form. Also, some solutions for pij in the last of our system of equations 
are purely hydrodynamic and these constitute the normal solutions [46]. 
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